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Abstract 

We reconsider the construction of the vector and axial-vector currents in Wilson Chiral Pertur- 
bation Theory (WChPT), the low-energy effective theory for lattice QCD with Wilson fermions. 
We discuss in detail the finite renormalization of the currents that has to be taken into account in 
order to properly match the currents. We explicitly show that imposing the chiral Ward identities 
on the currents does, in general, affect the axial-vector current at 0(a). As an application of our 
results we compute the pion decay constant to one loop in the two flavor theory. Our result differs 
from previously published ones. 
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I. INTRODUCTION 



Chiral perturbation theory (ChPT) is widely regarded to be an important tool for Lattice 
QCD. It provides analytic guidance for the chiral extrapolation of the lattice data obtained 
at quark masses heavier than in nature. Standard ChPT, as formulated in Refs. [ll, 0], is 
based on the symmetries and the particular symmetry breaking of continuum QCD. The 
generalization to Lattice QCD with Wilson fermions, taking into account the explicit chiral 
symmetry breaking of Wilson's fermion discretization [sl , was given in Ref . 01 . The resulting 
low-energy effective theory, often called Wilson chiral perturbation theory (WChPT), is a 
double expansion in the quark mass and the lattice spacing, the two parameters of explicit 
chiral symmetry breaking. 

Massless continuum QCD is invariant under various non-singlet chiral transformations. 
This invariance implies the existence of conserved currents (which are obtained by the 
Noether theorem) and various chiral Ward identities. ChPT is constructed in such a way 
that these Ward identities are correctly reproduced, order by order in the chiral expan- 
sion. And since conserved currents do not renormalize it is straightforward to maintain the 
normalization of the currents. 

The construction of WChPT is slightly more complicated compared to continuum ChPT. 
Due to the explicit breaking of chiral symmetry by the Wilson term there does not exist a 
conserved axial- vector current for vanishing bare quark mass. And even though a conserved 
vector current exists for degenerate quark masses, it is often not used in practice. The local, 
non-conserved vector current is employed instead, even though it requires the computation 
of a renormalization constant Zy. The renormalization constant Za is also needed for the 
axial- vector current. 

The explicit breaking of chiral symmetry and the need for renormalizing the currents 
raises the question how to construct the effective currents in Wilson ChPT. The "Noether 
link" does not hold anymore. Also the renormalization of the currents has to be taken into 
account for a proper matching of the effective theory to the fundamental lattice theory. 

Some results concerning the currents can be found in the literature [s], 0] , but they are 
in conflict. Reference [sl calculates the pion decay constant using the current obtained by 
the naive Noether procedure as the axial- vector current 0]. Reference introduces source 
terms for the currents as in continuum ChPT, and constructs the generating functional 
for correlation functions of the currents. The resulting axial-vector current contains an 
additional 0(a) contribution that is not present in the Noether current. Consequently, the 
resultant /,r contains an extra term and differs from that of Ref. [sl. 

Besides this discrepancy the issue of renormalization has not been properly taken into 
account in either work. No particular renormalization condition for the axial-vector current 
has been imposed as is necessary for a proper matching of the currents. It has been argued 
in Ref. [sl that the results for the currents derived there should hold for any choice of lattice 
operators which are correctly normalized in the continuum limit. However, the validity of 
this expectation has not been shown so far. 

In this paper we reconsider the construction and mapping of the vector and axial- vector 
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currents in WChPT to 0(a). ^ We proceed in two steps. First, we write down the most 
general expressions for the currents which are compatible with locality and the symmetries 
of the underlying lattice theory. With this procedure we reproduce the results of Ref. [gI 
(which are more fully justified in Ref. j^). In the second step we impose the chiral Ward 
identities as particular renormalization conditions for the currents. This choice, suggested 
m Refs. 0, lUl, [12], is widely used in practice. We find that this renormalization condition 



does have an impact at 0(a) on the axial- vector current. Consequently, our current differs 
from the ones in Refs. [sl, |6|. As an application of our results we finally compute the pion 
decay constant to one loop, including the 0(a) correction to the chiral logarithm, which also 
differs from the results in 0] . 

This paper is organized as follows. In section [TTl we first summarize some definitions of the 
lattice theory with two flavors of Wilson quarks, in particular the various vector and axial 
vector currents used in numerical simulations. This is followed by the Symanzik expansion 
of the currents close to the continuum limit. The currents in the Symanzik effective theory 
are then mapped to their counterparts in the chiral effective theory. Section IIIII discusses 
the renormalization of the vector and axial-vector currents in the lattice theory and how 
this is carried over to the effective theory. The results for the decay constant are given in 
section HVt followed by some concluding remarks in section |Vl Appendix |X] is devoted to 
an alternative but equivalent derivation of the currents based on the generating functional, 
while appendix [B] consists of details concerning the calculation of Za in the effective theory. 



II. CURRENTS IN WCHPT 

A. Definitions in the Lattice theory 

We consider Lattice QCD with Wilson fermions on a hypercubic lattice with lattice 
spacing a. For simplicity we study A'^^ = 2 quarks with equal quark mass. The fermion 
action is of the form 

Sf = S'w + CsW'S'clover , (l) 

where the first part denotes the standard Wilson action with bare quark mass ttiq. We 
also allow for a clover- leaf term with coefficient csw- The details of the gauge action are not 
important in the following so we leave it unspecified. 

It is common to use the local expressions for the vector and axial vector currents in 
numerical simulations, 

V^,Loci^) = ^(a;)7^T>(a:), (2) 
A';,Loci^) = ^{x)^,l,Ty{x). (3) 



Preliminary results have already been presented in Ref. [8| . Details have changed but the overall conclu- 
sions are unaltered. 
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The T"' are the hermitian SU(A^/) generators, normahzed according to tr (T°T'') = S""^ /2. In 
the case of Nj = 2, the one we are considering, this normahzation corresponds to = cr"/2, 
where a" are the usual Pauh matrices. 

For degenerate quark masses the fermion action (fl]) is invariant under S\]{Nf) flavor 
transformations. The associated conserved vector current differs from (121) and reads llOll 



^Con(^) 



1 f- 



+ ip{x + afi)U ^{x)il){x) + ii{x)Ul^{x)ip{x + afi) \ . 



(4) 



No conserved axial vector current exists due to the explicit chiral symmetry breaking by the 
Wilson term in S'w- 

In on-shell 0(a)-improved lattice theories with degenerate quarks one defines improved 
currents by adding terms involving lattice derivatives to the local currents 13 , 
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where 



^lmp(^) 



^^,Imp(-^) 



(1 + byam) 
= (1 + hxam) 

i){x)ia^yT''%l){x) 



^mVc(^)+cv^(V+ + V.)T;,(x) 



^:,Loc(a:) + CA7T(V^ 



and 



(5) 
(6) 

(7) 



The coefficients 6v,a and Cv,a — together with csw — can be non-perturbatively tuned such 
that the cut-off effects are of O(a^) instead of linear in a?' 

In order to correctly approach the continuum limit the non-conserved currents need to 
be renormalized. Thus one introduces 



' /i,ren,Loc 
/j,rcn,Loc 



^V,Loc ^/^^Loc 1 
^A,Loc ^)l,Loc ' 



V T - 
' /i,ren,imp 

/j,ren,Imp 



^V,Imp ^Imp 5 
■ ^A,Imp ^;x,Imp • 



(8) 
(9) 



In the following, we will often use 2^v,a generically, without specifying the underlying current. 

The Z-factors (which depend not only on the choice of currents but also on the action) 
can be fixed by imposing chiral Ward identities 10|, lUl, ll2|. We will come back to this 



important issue in section UTTl The Z-factor for the conserved vector current is 1, of course. 



B. The Symanzik effective theory 

According to Symanzik the lattice theory can be described by an effective continuum 
theory provided one is close to the continuum limit 15|, ll6| . This effective theory is defined 



^ Note that the mass m in eqs. ([5|) and ([6]) denotes the renormahzed mass containing the additive mass 
renormahzation proportional to 1/a and the renormaUzation factor: m = Zm{mo — Trier) /a. The factors 
of (1 + fcv.Aom) can also be considered to be part of the renormahzation factor, but it is notationally 
convenient here to include them in the bare currents. 
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by an effective action and effective operators, and botli are strongly restricted by tlie locality 
and the symmetries of the underlying lattice theory. The leading terms are, by construction, 
the familiar expressions of continuum QCD. Lattice artifacts appear as higher dimensional 
operators multi plie d by appropriate powers of the lattice spacing. The effective action, for 



example, reads [17 



5'sym = Sct + acsw J d x^{x)ia^yF^i,{x)^{x) + 0{a ). (10) 

The first term is just the continuum QCD action for two flavors with degenerate quark mass. 
The leading cut-off effects are described by a single correction term, a Pauli term containing 
the field strength tensor Ffj,y{x) multiplied by an unknown coefficient ('low-energy constant') 
csw Many more terms appear at O(a^) jlTI ]. 

The mapping of the bare local currents of of eqs. ([2]) and ([3]) into the Symanzik theory is 



13, 14 



^m"loc ^ ^;Sym,Loc = + ^V^^) (V;';,, + aCy^.T;,,,,) + 0{a^) , (11) 

^M,Loc ^ ^^,Sym,Loc = + ^A«"^) (^^,ct + "CA^^Pet) + 0{a') (12) 

^Loc ^ ^sWoc =^{l + hpam)P:, + 0{a) (13) 
r;.,Loc ^ ^;.,Sym,Loc = 4(1 + &Tam)T;,_,, + 0{a) (14) 
where the continuum bihnears take their usual forms 

^^ct = ^7m2^> , r;,,ct = ii^,.T^^ , Ai^, = i^^.^.ry , p^, = ^^.ry . (15) 

The mapping of operators between effective theories, 

OLat ^ Osyn., (16) 

is defined so that 

(OLat(x)0(^^^\ Al^\ y))s,^, = (Osy^(a:)0(^Sym^ ^Sym^ ^Sym^ ^^^^ (^7) 

Here 0{ilj^''\ip^''\ A^''\y) and 0(?/;^^™, -j/^^^™, A^^™, y) are arbitrary multi-local operators 
consisting of quark and gluon fields at positions y = yi, 2/2, ■ ■ ■ which all differ from x. The 
above vacuum expectation values are calculated with ^Lat on the left-hand-side and Ssym 
on the right- hand- side. We note that the equations of motion have been used in order to 
reduce the number of 0(a) terms on the right hand sides of f|T0l) - f|T2|) . This is legitimate as 
long as we consider the Symanzik theory as an on-shell effective theory, set up to reproduce 
physical quantities like masses, decay constants, etc.. 

We stress that the currents V^"sym,Loc ^^,sym,Loc result of matching the bare 

lattice currents into the Symanzik theory, and thus must include the renormalization con- 
stants Zyj^. These have perturbative expansions of the form Zyj^ = l + 0[g{ay], but do not 
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contain any contributions of 0(a), since all 0(a) terms are explicitly accounted for. They 
are completely determined in principle once one has specified the action and g{a), but are 
only known approximately in practice. This will not present a problem, since we ultimately 
will normalize the currents non-perturbatively and all dependence on Zy^ will cancel. The 
superscript "0", which indicates that this quantity is of O'th order in an expansion in a, 
distinguishes these Z-factors from the non-perturbatively determined renormalization con- 
stants Zv,A introduced before. The latter, as we will see, depend linearly on a. 

We can also map the renormalized currents of eqs. ([H]) and ([9]) into the Symanzik the- 
ory. For V^,rcn,Loc and A^^ren.Loc, the result will be the same as in (ITTi) and (|T2l) except 
for multiplication by overall factors of (the to-be-determined quantities) 2'v,loc and 2'a,loc, 
respectively. 

The results f lTTjl and f|T2l) also hold for improved currents of the form of eqs. ([5]) and ([6]), 
although the coefficients Zy ^ and cv,a will differ. It is important to keep in mind that these 
coefficients, as well as csw? are a priori unknown, their values depending on the details of 
the lattice theory. If all these parameters vanish we say that the lattice theory and current 
matrix elements are 0(a)-improved. For this to happen the parameters cx and bx of the 
lattice theory need to be tuned to appropriate nonzero values.^ 

The vector current in ffTTj) is the effective current into which the local lattice current is 
mapped. The conserved lattice vector current is mapped onto the most general conserved 
effective vector current. This current V'^sym.Con "^^^ t)e obtained by starting from (ITTl) (with 
an a priori different coefficient cy in front of the tensor term, and with = 1, 6v = 0) and 
then imposing current conservation. The current in (ITTl) is, however, already conserved, 

^M^M^Sym.Loc = , (18) 

which is a consequence of the particular structure of the 0(a) correction (total derivative of 
an antisymmetric tensor). Hence, both lattice currents are mapped onto the same form of 
effective current at this order in the Symanzik expansion. Violations of current conservation 
which imply a difference between the local and the conserved current are expected to show 
up at O(a^). 

We mention two special cases of the Symanzik effective theory at 0(a) that we will need 
in the next section for the mapping to the chiral effective theory. Let us first consider the 
Symanzik effective theory in which csw is non-zero but cv,a and 6v,a vanish, and Zy^ = 1. In 
this case the Symanzik currents coincide with the continuum QCD currents. These transform 
into each other under infinitesimal SU(2) axial rotations Sip = uj"'T°''-f5ip, Sip = u^'ip'y^T"': 

This leads to various chiral Ward identities, schematically written as (^^symC^Sym) = 
(SOsym), where Osym denotes some product of vector and axial- vector currents. The 
form of these Ward identities is as in continuum QCD, since the right hand side reads 
(^C^Sym) = (^C'ct)- The only difference is the appearance of an extra term proportional 



^ Explicitly, cy = -{Z^/ Z^)cy, ca = -{Zp/ Z'^)ca, by = -by, and 6a = — &A- 
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to acsw in the variation of the effective action, caused by the Pauh term in 5'sym- Note, 
however, that these simple QCD-like Ward identities only hold at 0(a), and are violated at 
O(a^) by the terms of this order in the effective action and effective currents. 

The second special case is obtained by setting csw = 0, 6v,a = and ^va — 1' 
which case 0(a) corrections stem entirely from the cv,a terms in the effective currents. This 
implies, for example, that the correlation function of two axial-vector currents reads 

(^;:,Sym(^)^",Sym(y)) = (^^,ct(^)^:,ct(y)) 

+acA{A;^A^)d,P:tiy) + d,P:,{x)Al^,,{y)) . (20) 

Here the expectation values are defined with Boltzmann factor exp{—Sct) only. Hence, the 
0(a) correction is given by the correlation function between the axial- vector current and 
the derivative of the pseudoscalar density. Again, this result will be violated as soon as one 
includes the corrections of O(a^). 



C. Matching to ChPT 

The appropriate chiral effective theory is obtained by writing down the most general chiral 
effective Lagrangian and effective currents which are compatible with the symmetries of the 
underlying Symanzik theory. A standard spurion analysis is employed in order to properly 
incorporate explicit symmetry breaking terms. For example, the Symanzik effective action is 
invariant under the chiral symmetry group G, parity P and charge conjugation C, provided 
both the mass and the coefficient acsw are promoted to space-time dependent external fields 
M and A, which are postulated to transform according to 

M LMR\ M Aft, M ^ , . . 

A^LAR\ A^A\ A^A^. ^ ' 

The "physical" values are obtained by setting M ^ m and A — > acsw- In an intermediate 
step, however, the spurion fields M and A are used together with the standard field 

S(x) = exp (^Ti-{x)T'^^ (22) 

in order to write down the most general scalar that is invariant under G, P and C. This 
has been done in Refs. [s], [isj , and part of the result reads (in Euclidean space-time) 

>Cchirai = y (9^S(9^S)t) - ^2Bm{J:^ + S) 

+L4525m(a^S(a^S)t)(St + S) 
+W^45acsw(5^S(9^S)t)(St + S) 

-We82Bmacsw{^ + + ■■■ ■ (23) 
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Here the angled brackets denote traces in flavor space. ^ The lattice spacing appears in the 
combination 

a = 2Woa, (24) 

which is of dimension two.^ We have dropped a number of terms of 0{p^) and O(a^), which 
we will not need in the following. Note that we have absorbed the term pWoOCswi^'^ + ^)/2 
in the deflnition of the mass, so m in fl23l) corresponds already to the so-called shifted mass 

0,0. 

Taking the naive continuum limit a ^ we recover the correct terms of the continuum 
chiral Lagrangian with the familiar low-energy coefficients /, B and L45 = L4 + L5/2 of 
continuum ChPT [l|, [2|- Wq and = W4 + W5/2 are low-energy constants associated with 
breaking terms due to the nonzero lattice spacing (sf.^ 

We now apply the same procedure to derive expressions for the effective operators. The 
currents in the Symanzik effective theory are given in eqs. fill I) and fll2l) — forms which, as 
noted above, hold for all the lattice currents of interest. To simplify the following discussion, 
we will map the Symanzik currents without overall Z-factors and bx corrections into the chiral 
effective theory. These factors can be added at the end. Thus we consider the mappings 

i;^,, + acv5.T;,,et — ^^efr , (25) 
A;^,, + a-CAd.P:, — > . (26) 

The leading contributions are just the continuum expressions for the currents, 

^^ct = ^«7mT>« + ^lI.T^^Pl , (27) 
^M,ct = i^R7,Tyn - ^L^.T^L, (28) 

where we decomposed the currents into chiral fields. The vector current is just the sum of 
the right and left handed current, the axial- vector current is given by the difference. 

The 0(a) corrections couple fields with opposite chirality, and the currents in the 
Symanzik theory do not transform under chiral rotations as continuum vector and axial- 
vector currents. However, the continuum transformation behavior can be enforced by pro- 
moting the coefficients cy,ca to spurion fields Cy and Ca with non-trivial transformation 
behavior under G, P and C: 

Cx ^ LCxR\ Cx ^ C],, Cx^Cl, X = V,A. (29) 

Note that these transformation laws are identical to the ones of the other 0(a) spurion A, 
cf. dUl). 



^ In the last section wc used the same notation for correlation functions. The context usually tells unam- 
biguously what is meant by (. . .). 
^ Wo, a LEG that enters the chiral Lagrangian at 0(a), is of dimension three [5|- 

® Our notation for the low-energy coefficients differs slightly compared to the notation in other references, 
since our m is already the shifted mass: our W45CSW and Wgscsw correspond to the combinations W45— L45 
and Wes — 2L68 in Ref. These combinations are abbreviated to W and W in Ref. 0. 
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It is now easily checked that the 0(a) corrections written in the form 



+ {d,i^j,)ta^,ClT^ijL + i^nta^.T'^Cld^ijL, (30) 

+ {d,i^j,h,C\T^ijL + i^nl^T'^Cl.d.^L, (31) 

transform as the continuum currents. Setting the spurion field to its "physical" value, 
Cy — > acy and Ca clca, one recovers the desired Symanzik currents. 

The currents in the chiral effective theory are now obtained by writing down the most 
general vector and axial- vector current, build of the chiral field S, its derivatives and the 
spurion fields. The spurious necessary for the construction of the vector current are M and 
A, the ones we already encountered in the construction of the effective action, and Cy- For 
the axial-vector Ca must be used instead. 

In order to write down the currents we first recall that our effective theory has to reproduce 
continuum ChPT if we send a ^ 0. This requirement implies that the continuum part of 
the currents are just the expressions given by Gasser and Leutwyler. At leading order these 
read 
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V;,LO = y (T'^(St9^S + S9,,St)) , (32) 



^M,LO = y (T'^(St9^S - E9^St)) . (33) 

Obviously these expressions transform as vector and axial-vector currents under G, P and 
C. Moreover, these currents are properly normalized in order to satisfy the current algebra. 

In order to construct the leading 0(a) corrections we need at least one power of either A 
or Cx and one partial derivative of S. It is easily checked that the following terms transform 
as desired: 



VI) y;Lo(s^^ + sAt), 

V2) V^l^oi^^Cv + ^Cl), 

V3) {T'^id^J^^A - A9^St + d^EA^ - A^d^E)), 



(34) 



for the vector current, and 



Al) A;i,Lo(SM + SAt), 

A2) A;i,lo(S^C^ + SC1), 

A3) {T^id^E^A + Ad^E^ - d^EA^ + A%E)), 

A4) {T'^{d,j:^CA + CAd^^^-d,j:c\ + cl,d,j:)), 



(35) 



for the axial- vector current. Setting the external fields to their final value, Cx o-cx, we 
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obtain the following expressions for the currents in the effective theory: 



ya 



(W^yiCsw + W^y2Cv)(S + St)^ , 

a(W^AiCsw + W^A2CA)(S + St) 
+4a(iy^3Csw + WAiCp:)d^{T\T. - T))) . 
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(36) 



(37) 



The coefficients Wx are unknown low-energy constants (LECs). In order to make Wxi-, Wx2 
dimensionless we included the factor 4//^ . Note that the second line in ( l37l) is proportional 
to the continuum pseudoscalar density, P^^ = pB{T"-(T, — T))) /2. 

The number of unknown LECs in the currents can be reduced using the freedom of field 
redefinition j^, 19]. Explicitly, performing S ^ S + (5S with 
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<5E = — cswAW^(s2-l), 
we obtain the same effective Lagrangian and currents with the transformed coefficients 



(38) 



Wvi 

Wai 



1^45 + , 

Wvi + Aiy , 
Wai + Aiy , 



Waz ^ Wa^-^W. 



(39) 



Therefore, depending on the particular choice for AVT we may make one LEG vanish. In 
the following we choose this to be Waz in the expression for the axial- vector current. 

So far the expressions in (l36ll and ( 1371) are the most general currents compatible with the 
symmetries that have the correct continuum limit. In order to match the currents properly 
we have to impose the constraints that these currents have to obey at 0(a), for instance 
current conservation for the vector current. This will relate some of the LECs Wx to those 
in the effective action. 

In order to discuss this we first derive the EOM corresponding to the Lagrangian in eq. 
f l23|) without the continuum NLO term proportional to L45 and without the 0(am) correction 
proportional to 1^68-'^ Using the shorthand notation x = '^Bm, p = 2W45WQacsw aiid 



R 



4 



j2 



(40) 



the leading order EOM reads 



S(9^9^St) - (9^9^S)St (1 + pi?) 



2A 



S - St 
pdy,R. 



ix + pT) 



(41) 



The parameter A is the Lagrange multiplier associated with the constraint det S = 1. Setting 
the lattice spacing to zero (i.e. p = 0) eq. (HTi) reproduces the EOM in continuum ChPT Q]. 



^ These can and need to be included if the NLO terms of 0{p^m) and O(p^ma) are included in the current. 
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Using (jHj) we find the condition 



d^V;^ = ^ Wvi = W45 and Wv2 = ^ ■ (42) 
Therefore, the conserved vector current in the chiral effective theory is given by 

^M^fT = V^^^o (1 + -^aW^45Csw(S + St)^ . (43) 



This expression agrees with the one in Ref. [6|] obtained from the generating functional 
Note that at this order this current coincides with the Noether current associated with 
vector transformations. This will, however, no longer be true at higher order in the chiral 
expansion.^ 

Note that there is no term proportional to cy in ( H3i) . which is a consequence of the d^T'^y 
structure of the 0(a) correction. We need three partial derivatives in order to construct 
such a term in the chiral effective theory. Hence, this correction is of O(ap^), which is of 
higher order than we consider here. 

In order to obtain the proper result for the axial-vector current we have to make sure 
that the properties (fTOll and ( l20l) of the underlying theory are correctly reproduced. The 
result ( l43l) for the vector current together with ( fT9|) (for cy = ca = 0) immediately implies 
Wai = W45. On the other hand, setting csw = and demanding (120|) leads to Wa2 = 0. So 
the final result for the axial vector current reads 

^M,eff = ^Uo (1 + -^aW^45Csw(S + St)^ + AaWACAd,{T^iJ: - St)) , (44) 



where we abbreviated Waa = Wa- This expression also agrees with the result in Ref. [6|.^ 

So far we have discussed the currents at leading order in the chiral expansion including 
the first correction of 0(a). Higher order contributions to the currents can be derived in the 
same fashion. The terms without factors of the lattice spacing are the familiar contributions 
of continuum ChPT, 

= y (T"(St9^S + SS^St)) (J^2BmL,,{J: + St)^ , (45) 



(46) 



In order to construct the first subleading 0(a) corrections we have to form vector and axial- 
vector currents with one power of either A or Cx, and either three derivatives (corresponding 
to the 0(ap^ contributions), or one derivative and one power of the mass spurion field M 



At higher order in the chiral expansion the currents receive contributions which are not present in the 
Noether current. The vector current contribution proportional to Lg Q, which captures the dominant 
contribution of the pion form factor, is an example for this. We thank J. Bijnens for pointing this out to 
us. 

The result in Ref. Q is obtained by setting Wio = 2Waca and W = 2W45Csw- 
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(the 0(am) contributions). Simple examples for such terms are products of the leading 0(a) 
terms in (!34l) - (l35|) with either (EM''' + MS"!") or {d^'E{d^T,y). An example is the correction 

6A;^,^ = ^(T'^(St9^S - S9^St)) (XcAam(S + St)^) , (47) 

with X being the low-energy constant associated with this correction. 

There are more terms possible and it is a straightforward but tedious exercise to find all 
possible terms contributing to the currents. For the rest of this paper, however, we do not 
need these higher-order terms. 

We conclude this section by discussing the pseudoscalar density, which we need later on 
for the computation of the PCAC mass. In the Symanzik effective theory the pseudoscalar 
density is given by P^,". The corresponding expression in the chiral effective theory is con- 
structed in the same way as the currents. Since the density in the Symanzik theory has no 
0(a) term, we only have the spurion A for the construction of the density, and we find^° 

= ^{T''^ - S)) (l + ±aW,sCsw{^ + S^)^ . (48) 

As already mentioned, our expressions for the effective currents agree with the results in 
Ref. [sl, although the two calculations use different methods. In the approach of Ref. j^, 
which follows the procedure used by Gasser and Leutwyler in continuum ChPT [ij, sources 
for the currents are introduced and the generating functional is constructed. The currents 
are then obtained from the generating functional by taking derivatives with respect to the 
sources. The analysis is complicated, however, by the presence in the Symanzik theory of 
0(a) violations of the local gauge invariance used to restrict the mapping to the chiral effec- 
tive theory. This complication, missed in Ref. j^, was noted in Ref. 9|, and an appropriate 
extension outlined. The outcome of this extension was that the form ol there^ute ot Ref. 
still hold after a redefinition of the LECs. Since a systematic and complete treatment of 
the impact of the 0(a) effects in the generating functional method has not been presented, 
however, we provide such a treatment in appendix |Al In particular, we show that the same 
results ( 143|) and (144|) are obtained using this method. 



III. RENORMALIZATION 

We now return to the issue of the normalization of the currents. The results of the 
previous section allow us to determine the form that a given lattice current will have when 
mapped into WChPT. For example, putting back the overall factors, the renormalized local 
vector and axial- vector currents of eqs. ([8]) and iQ map as 

^M,ren,Loc ^ ( 1 + &Vam) 1^% , (49) 

Z — 

^/.,ren,Loc ^ -^{1 + bAam)A^^^^ff , (50) 



This agrees with the resuh of Ref. 3i with Wqs related to W of that reference hy W = WqsCsw- 
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with F^gfj and A^ eS given by eqs. (H3|) and (jS]), respectively. In the lattice theory, the 
renormalization factors are determined by imposing particular conditions at nonzero lattice 
spacing. Hence, in order to properly match the effective currents we should impose the same 
conditions in the effective theory. 



Renormalization of the lattice currents 



Since the conserved vector current does not need to be renormalized it can be used to 



normalize the local current and to define ^v,Loc (af the lattice level) according to 20 



21 



^ (/iv;%nK) 



(51) 



Here i and / denote arbitrary initial and final states, though it is convenient to choose 
zero-momentum pseudoscalar states. Note that the dependence of ^v,^ on the particular 



states can be sizable, in particular if the theory is not 0(a)-improved 20 



An alternative definition for Zv,loc, which does not use the conserved current, is given 



by [20 



21 



22 



^V,LocK(p)|<Lock^(p)) 



(52) 



The matrix element on the left- hand- side can be obtained in the usual way by calculating 
the ratio of two correlation functions, where pseudoscalar sources are used to project onto 
the pion states. Although we have specified initial and final pion states, the renormalization 
factor still depends on the momenta of the pions 20|] . For massive pions one usually chooses 
zero spatial momentum, p = 0. 

The two renormalization conditions ( 15T1) and (l52l) do not specify the renormalization 
condition completely, since the matrix elements still depend on the quark mass. Theoreti- 
cally preferable are mass independent renormalization schemes, in which the renormalization 
condition is imposed at zero quark mass. This implies some technical difficulties, because 
numerical simulations cannot be performed directly for vanishing quark masses. One way 
to circumvent this technical limitation is to calculate the Z-factors for various small quark 
masses and extrapolate the results to the massless limit. This may introduce some extrap- 
olation error, but in principle is a viable procedure. 

A different strategy to define and compute the renormalization factors makes use of 
Schrodinger functional boundary conditions [23|. With this setup it is possible to compute 
■^v,Loc and other renormalization factors for a vanishing partially conserved axial-vector 
current (PCAC) quark mass. Even though the renormalization conditions are imposed at 
vanishing quark mass, they now depend on details of the Schrodinger functional setup, e.g. 
the size and geometry of the finite volume. This dependence can, in principle, be removed 
by extrapolating to infinite volume. 

Being able to work at vanishing quark mass has another adva ntag e: chiral Ward identities 
involving the axial- vector current simplify significantly. In Ref. 24], for example, the Ward 
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identity 



abc I 



dR 



(53) 



is imposed. This identity is the Euchdean analogue of the current algebra relation stating 
that the commutator of two axial- vector currents is equal to the vector current 25|]. The 
region R is chosen to be the space-time volume between two hyper-planes at = i/o ± 
t. The equation for z/ = between pseudoscalar states has been used to determine the 
renormalization factor Z^. For more details see Ref. 24 . 



B. Renormalization of the effective currents 

Having chosen particular renormalization conditions for the lattice currents, we have to 
impose the same conditions in the chiral effective theory. 

Some conditions are harder to implement than others. For example, matrix elements 
between the vacuum and a vector meson state in eq. (ISTl) are not easily accessible in standard 
mesonic ChPT, since the vector meson is not a degree of freedom in the chiral effective 



theory. Conditions involving quark states (the so-called "RIMOM" scheme [26|) are also 
out of reach. In practice, only conditions involving pseudoscalar states can be treated in the 
chiral effective theory. 

In the following we carry out the matching for one particular class of renormalization 
conditions. For the vector current we assume that either the condition flHTl) or (jS^D is 
imposed with single p ion st cit G s cit zero spatial momentum and at vanishing bare PCAC 
mass. For the axial-vector current we assume that condition (1531) is employed to fix Za- 
We impose these conditions in infinite volume. Finite volume can also be considered, but it 
does not make a difference at the order in the chiral expansion to which we work. 

As a first step we have to calculate the PCAC mass and set it to zero. The PCAC mass 
is defined by 

m.cAc= ^^,\pa^\,.^ ■ (54) 

Expanding the S fields in (jHj) and (H8l) . keeping only one power of tt", the ratio of expectation 
values on the right hand side is easily obtained^^ 



Ml 

2B y ' p 



mpcAc = 7T^ ( 1 + 7^a[2(iy45 - Vl^68)csw + WaCa] ) , (55) 



In practice, to calculate mpcAC requires knowledge of the renormalization constants of the lattice axial- 
vector current and pseudoscalar density, which can. as the present work shows, introduce additional 0(a) 
corrections. These do not, however, change the key result being derived here, namely that mpcAC oc m 
up to O(a^) corrections. 
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in agreement with Ref. [6|]. The PCAC quark mass is proportional to the pion mass, which 
is given by [gI 

= 2Bm (^1 + ^a{2We8 - W^45)csw) • (56) 

Recall that m denotes the shifted mass including the leading 0(a) shift. 

At higher order in the chiral expansion the right hand side of fl56|) contains an additional 
correction of O(a^) 0]. A vanishing pion mass therefore corresponds to m = O(a^). Since 
we ignore O(a^) corrections, we conclude that a vanishing PCAC quark mass is equivalent 
to m = 0, and in the following we assume this condition for the shifted mass.^^ 

The next step is the determination of Zv,loc using either (|5T1) or fl52|) . Both conditions 
are easily calculated using the chiral effective theory "image" of the local lattice current, 
eq. fl49l) with m = 0, and the expression fl43p . In both cases we find^^ 



^V,Loc 



Z'y . (57) 



Even though the result is the same for both renormalization conditions, the way it arises 
differs in the two cases. The result for the first condition is obviously trivial since both the 
local and the conserved effective vector currents have the same form in the effective theory 
[cf. ( 143|) ]. differing only by overall factors. In (1521) . on the other hand, the two pion states 
contribute a wave function renormalization factor Z^, which reads^^ 

16 

= 1 - j^aWi^csw, (58) 

and which exactly cancels the 0(a) correction coming from the current. 

Let us consider higher order corrections to the result ( l57j) . Expanding the S fields in the 
vector current (H3l) the first correction terms contain four pion fields. Two of these need to 
be contracted and form a loop. Hence the result will be proportional to M^lnM^, which 
vanishes for mpcAc = 0. Higher order analytic parts, on the other hand, will contribute 
corrections of O(a^), which we do not consider here. 

The implication of eq. (!57|) is that, after one or other of the renormalization conditions 
has been enforced, the renormalized local vector current maps simply into V^eff' since the 
prefactors on the right-hand-side of eq. (149!) cancel. The same holds true for any lattice 
vector current which is renormalized in this way. 



The O(a^) correction implies a non-trivial phase structure of the theory with two qualitatively different 
scenarios [J| . One of these is characterized by a first order phase transition and the pion mass is non-zero 
for all values of to. For to = 0, however, the pion mass assumes its minimal value ^-^^^ = O(a^). Since 
here we ignore the O(a^) corrections we also find for this scenario that a vanishing PCAC mass is given 
by a vanishing shifted mass, at least to the order we are working. 

It is a simple matter to restore the mass dependence, in which case ([57| would read ^v,Loc = Zy{l+byam). 
Note that there are no terms proportional to the quark mass in the chiral Lagrangian (pS)) , since we have 
set to to zero. 
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We now proceed to the local axial- vector current, whose normalization ^a,Loc is to be fixed 
by imposing Ward identit y (l53ll . using external one-pion states having definite (nonzero) 
momenta. Following Ref. [2J] we choose the region R to be the space-time volume between 
two hyper-planes a.t xq = yo ± t with some finite time separation t. The equation for z/ = 



can be brought into the form 24 



dxe'^^^e'^^'^in^p) \ [Ag,_(yo + t,x)- A^^^^/o - t, x)]Aireniy)\^m 

= 2te^'^Ap)\Vo:rcniy)\^m , (59) 

The matrix element on the right hand side of this equation is essentially the one in (!52|) that 
we used to fix Zv,loc- Imposing (!59|) thus provides a condition for Za,loc: simply compute 
the two sides in the effective theory and set them equal. The calculation is straightforward 
but rather technical. For this reason we defer the details of the calculation to appendix [Bl 
The final result is (recall that m = so there is no 5a term) 



JA,Loc 



4d 

- l-^{W^5CsW + WACA)zA{t), (60) 

ZA{t) = l-cosh[t{\p\-\q\)]eM-\t\\p-q\]. (61) 

Since this result is determined by a ratio of physical correlation functions in which the axial 
currents are separated, it must depend on the "physical" combination of LECs W45CSW + 
WaCa • This point is also explained at the end of appendix \^ and provides a check of 
our result. We emphasize that fl60|) is the complete result to 0(a), there are no corrections 
of 0(am) since we work at zero quark mass. The next correction to (I6OI) is of O(a^), which 
we do not consider here. 

In contrast to the vector current we do find a non- vanishing correction of 0(a). That this 
correction depends on the separation between the axial currents, t, and upon the external 
state, is as expected. The t dependence is proportional to an integral (a sum on the lattice) 
of the divergence of the axial current, which does not vanish at 0(a), even when m = 0. The 
dependence of the external state is a generic feature of 0(a) corrections in an unimproved 
theory. 

Using eqs. ( |50l) and (|60l) we find that the renormalized local axial-vector current maps 
into WChPT as 



AO. ^ AO, 

/i,ren,Loc /x,ren,Loc,cff 



4d _ 

1 - -p (W^45CSW + Wj(Cp^ZA{t) 



^M,efr • (62) 



This is the main new result of this paper. We see that the a dependence of ^^ eff' derived 
using symmetries and given in eq. (jUj) , is supplemented by an additional discretization error 
resulting from the application of the normalization condition at non- zero a. 

Note that the quantity Z\ appears only in the combination Za,loc/-^a ^qs. ( l50l) and 
f l60|l . and the individual value of Z\ is not necessary. This factor is needed only when express- 
ing the bare lattice current in the intermediate Symanzik theory. In fact, the combination 
^A,Loc/^A5 also the analogue for the vector current, Zv,loc/^V' may be interpreted as 
renormalization constants ^A,eff and 2'v,eff in the chiral effective theory [8|. 
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We also remark that the form of ^^^rcn,Loc,cfr apphes to any lattice current — local or 
improved. These cases only differ in the values of the LECs. Thus in the following we will 
drop the subscript "Loc" on the renormalized current in the chiral effective theory. 

We close this section by investigating the dependence of ZA^t) on t, p and q. There turn 
out to be three distinct cases (recalling that p, 0): 

• 'p = q. This is the simplest case to implement practically, and leads to ZA_it) = 0. Thus 
it turns out that, in this case, there are no additional 0{a) terms introduced by the 
current normalization. 

• p parallel to q. Then, for \t\ ^ l/\p — q\, the product of cosh and exponential becomes 
1/2, and so za 1/2. 

• All other non- vanishing p and q. Here, for \t\ ^ ^/\p — q\, the exponential overwhelms 
the cosh and za I- 

We stress, however, that in both the second and third cases za depends on t for non- 
asymptotic values of t. 



IV. APPLICATION: PION DECAY CONSTANT 

At this point we have completed the construction and matching of the effective currents. 
Now we can proceed and compute correlation functions involving these currents. As a simple 
but important example we calculate the pion decay constant. Given the presence of the ZA{t) 
contribution, our result differs, in general, from previously published ones. 



A. Decay constant at tree level 

Expanding the renormalized axial- vector current in powers of the pion fields we obtain 
to leading order the expression 

^;:,ren,efr = = (l + {W,,Csw + WaCa) [2 - ^A(t)]) • (63) 

This depends, as required, only on the "physical" combination W45CSW + WaCa of LECs. 
We can immediately read off the following expression for the leading order decay constant: 

fn,LO = / (^1 + -^a {W^5Csw + WaCa) [2 - ZA{t)]^ ■ (64) 

We see that the 0(a) corrections depend not only on the form of the underlying lattice 
action and currents (through the LECs) but also on the choice of renormalization condition 
[through ZA(t)]. The ZA(t) term always reduces the size of these corrections, although by 
no more than a factor of two (which happens in the third case discussed at the end of the 
previous section). The decay constant is free of 0(a) corrections only if both the action and 



17 



the axial-vector current are improved, i.e. for csw = ca = 0, in accordance with what we 
know from the Symanzik effective theory. 

We want to comment on the origin of discrepancies between (164|) and previously published 
results for the decay constant. Reference [sf finds^^ 

/So = / (l + ^aW,,csv?j , (65) 

There is no correction proportional to WaCa since their calculation used the Noether current 
as the axial- vector current [3]. This misses the 0(a) correction in f|T2|) . and cannot be correct 
since the result fl65l) is 0(a) improved if only the action (and not the current) is improved. 

Reference ^ finds the same form as flMl) except without the ZA{t) contribution. These 
authors assumed that a non-perturbative renormalization condition had been applied, but 
did not include the impact of applying the condition at non- vanishing lattice spacing. 



B. Decay constant to one loop 

It is straightforward to compute the leading correction to the tree-level result ( 164|) . Ex- 
panding the axial-vector current in eq. (1371) to higher powers of the pion fields one obtains 
the one-loop contributions to /^r. The integrals that appear can be regularized using dimen- 
sional regularization. The counterterms for the divergences are provided by the tree-level 
contribution of the NLO terms in the axial- vector current, c.f. ( H6|) and (147|) . Even though 
we have not explicitly given all possible NLO corrections of 0(am,ap^), it is easy to con- 
vince oneself that all contributing tree-level terms with one partial derivative are of the 
form ifdfj,Ti°' ■ am. Expressing m by the tree-level pion mass according to (l56l) . we obtain 
the counterterm 



(66) 



For simplicity we have absorbed the coefficients cgw and ca, and the contribution propor- 
tional to ZAif), in the LEG VFas, since in practice these and the W coefficients are difficult 
to disentangle. The additional factor 1//^ is introduced for convenience, since it leads to a 
dimensionless coefficient Was- 

The final one-loop result for /^r is then given by 



1 + -^W^A2 



M2 In ^ + — 



L45 + J^^- 



A3 



(67) 



The coefficients are Wai = 4(W45Csw + W^aca)[2 — ZA{t)] (as above) and Wa2 = 4(W45Csw + 
WaCa)[1 — ZA{t)]- Note that both coefficients depend on the physical combination of LECs, 
as expected. 



Ref. Q quotes explicitly the 3- flavor result, which we changed to the corresponding 2-flavor result to make 
the comparison. 
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A couple of comments concerning (1671) are in order. In the continuum limit we reproduce 
the familiar result for /^r from continuum ChPT for Nf = 2 Away from the continuum 
limit the result is modified not only by terms of 0(a), but also by contributions of O(aM^) 
and 0{aMl\nMl). 

The coefficients PVai and are in general not independent. For example, the first case 
discussed at the end of the last section has PVai = 2PVa2 (for asymptotically large t values). 
For the third case we even find VFa2 = 0, so the coefficient of the chiral logarithm is free of 
0(a) artifacts. Since in this case also Wai assumes its minimal value this is the theoretically 
preferred renormalization condition. 

Except for the special case with = the coefficient of the chiral logarithm receives 
an 0(a) correction in the form of the factor [1 + aW a2 / P]- Consequently, the coefficient of 
the chiral logarithm is, in contrast to continuum ChPT, not a universal coefficient depending 
on / (and Nj) only, but on the (non-universal) lattice artifacts too. This fact has previously 



been stressed in Ref. 27 



Note that the combination L45 appears effectively in the lattice spacing dependent com- 
bination L'§^{a) = L/^^ + oWazI P ■ Determinations of L45 based on simulations at one lattice 
spacing only are potentially dangerous because the size of the contribution aWA-i/ P is a 
priori unknown. 



V. CONCLUDING REMARKS 

In this paper we have reconsidered the construction of the vector and axial vector currents 
in WChPT. Due to the explicit chiral symmetry breaking in Lattice QCD with Wilson 
fermions two aspects need to be taken into account which are not present in continuum 
ChPT. 

Firstly, the local lattice currents are not conserved, and in general they do not map onto 
the conserved currents in WChPT. In particular, the WChPT currents are not the Noether 
currents associated with the chiral symmetries, not even at leading order in the chiral expan- 
sion. The reason is that the currents in the Symanzik theory have 0(a) corrections which 
are not related to the effective action. 

Secondly, the matching of the currents needs to take into account the finite renormaliza- 
tion of the local lattice currents. In order to properly match the currents the same renor- 
malization conditions that have been employed for the lattice currents must be imposed on 
the WChPT currents. 

Depending on the particular choice for the renormalization conditions the expressions for 
the renormalized currents differ by terms of 0(a). As a result, the WChPT predictions for 
matrix elements of the currents are different as well. Consequently, a result for an observable 
like /tt should make reference to the renormalization condition one has adopted. This has 
to be so at some level, since the lattice data differs too depending on the condition one has 
chosen. What we find is that the dependence enters at 0(a). 

At a technical level, our result does not change the number of low-energy coefficients that 
enter into WChPT predictions at 0(a). In particular, as stressed in Ref. j6|, there are only 
two combinations of the LECs that can enter into physical quantities, allowing their 0(a) 
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corrections to be related. What changes is the nature of these relations, which now depend 
to the choice of normalization condition. 

The considerations of this paper apply more generally. Another example from WChPT 
is the ratio of the renormalization factors for pseudoscalar and scalar densities. This must 
be determined non-perturbatively, e.g. by enforcing Ward identities, and will presumably 
receive 0(a) corrections similar to those for the axial current, although we have not worked 
out the details. The combination mS is, however, protected by exact lattice WTIs. 

More generally still, our results emphasize the (perhaps rather obvious) point that non- 
perturbative renormalization conditions generically introduce additional discretization er- 
rors. These must be (and in practice usually are being) accounted for when extrapolating 
to the continuum limit. 
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APPENDIX A: THE GENERATING FUNCTIONAL METHOD 

In section Hi CI we constructed the effective currents in a rather direct way. We first wrote 
down the most general currents that transform as vector and axial-vector currents. In a 
second step we imposed appropriate Ward identities that these currents must obey, and this 
led to various constraints on the LECs in these currents. 

We already mentioned at the end of section III CI that this is not the way one proceeds in 
continuum ChPT [l|, [2|. Instead, there one sets up the generating functional for correlation 
functions involving the currents (and the scalar and pseudoscalar density), and matches this 
to the analogous generating functional in the chiral effective theory. In this way the currents 
are obtained by functional derivatives with respect to the sources. 

A crucial link in the matching is the invariance under local chiral transformations, which 
plays the role of a gauge symmetry. The local invariance implies that the sources enter the 
chiral Lagrangian in a very restricted way, namely in form of a gauge-covariant derivative. 

One may ask if the same procedure is also possible once we include the lattice spacing 
corrections. In order to answer this question we derive in detail the generating functional 
for the Symanzik effective theory and match it to the one in WChPT. We will see that we 
indeed obtain the same currents as the ones in P3l) and (jH]). However, we will also see that 
it is more complicated to maintain the invariance under local chiral transformations. 

The "generating functional method" has been used before in WChPT in Ref . jol . However, 
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as discussed in Ref. [9|, the generating functional of Ref. [6|] violates the invariance under 
local chiral transformations at nonzero lattice spacing; only the leading (continuum) part 
respects this symmetry. Consequently, the construction of the axial-vector current requires 
some care: part of the 0(a) corrections have to be mapped separately into the effective 
theory The final result found in for the currents turns out, however, to have the form 
given in Ref. j^, as we derive here in a systematic way. 



1. Generating functional in continuum QCD and ChPT 



In order to prepare our discussion we first give a brief review of the construction in 
continuum ChPT, mainly to introduce our notation. We start by defining a (Euclidean) 
QCD Lagrangian including a part that includes sources for the currents and densities, 

^ = ^QCB — ^-^^Source , (Al) 

where £qcd is the usual massless QCD Lagrangian, while the second part 

'^^Source = V^7m ["fii^) + 75«m(^)] + [s(x) + 75p(a;)] (A2) 

contains sources for the vector and axial-vector currents (w^, a^) and for the scalar and 
pseudoscalar densities {s,p). These sources are matrix valued fields, given by 

v,{x) = vl{x)T\ a,{x) = al{x)T\ 
s{x) = s''{x)T\ p{x) =p"(x)T", (A3) 

where are the hermitian SU(A'^^) generators, normalized according to tr(T'*T^) = (5"^/2. 
We are interested in is Nf = 2, for which T" = (J°'/2, with a" the usual Pauli matrices. 

After integrating over space time we obtain the action in the presence of the sources, 
S = J d^xC{x). By taking functional derivatives with respect to the sources we obtain the 
vector and axial vector current and the scalar and pseudoscalar densities: 

= i^ixh.T'^i^ix) = K«(x) , (A4) 



Sv^{x) 

'S^) = ^(^)^M^5T>(x) = A;{x) , (A5) 
= i^{x)Ty{x) = S^{x) , (A6) 



Ss^ix) 

. 6S 
f 

6p"-{x) 



^(x)75TXx) = P"(x) . (A7) 



The key observation is that the Lagrangian (lAip is invariant under local SU{Nf) x SU{Nf) l 
transformations, which act on the fermion fields according to 

ip{x) ijj'ix) = R{x)P+ip{x) + L{x)P^ip{x), 

i^{x) ^'{x) = ^{x)P+L'^{x) + ^{x)P_R^{x) . (A8) 
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The projectors are defined in the usual way, P± = (1 ± 75)/2, and project onto fields with 
definite chirality, 



tpR 



IpL 



Crucial for the local invariance is the non-trivial transformation of the source fields 



-^11 
P 
P 



s' + p' = L{s+p)R'^ , 
s — p = R{s — p)L^ . 



(A9) 



(AlO) 
(All) 
(A12) 
(A13) 



(For notational simplicity we drop from now on the argument x). The invariance is more 
easily seen if we express the Lagrangian in terms of left- and right-handed fields, 



Source 



^ij7M K + «m] '>Pr + V^l7m K - «/^] + V^L + P] ipR + '4^r[s- p] 'ipL ,(A14) 



which obviously is invariant under global transformations. Note that global transformations 
leave £qcd and £sourco independently invariant. Under local transformations the derivative 
part of £qcd produces extra terms which are cancelled by the derivative terms in ( lAlOl) and 
(lAllI) . The reason why this cancellation works is related to the fact that the currents are 
the conserved Noether currents associated with chiral transformations and hence stem from 
local variations of the derivative term in the Lagrangian. 

The Lagrangian C can be conveniently rewritten in a form that makes the local invariance 
more transparent. Firstly, the source term in (]A14I) suggests to introduce sources for right- 
and left-handed currents. 



(A15) 



which transform according to 



r; = Rr^R^ + iRd^R^ 
I' = LLL^ + iLd^L^ . 



Vector and axial-vector currents are then obtained by 



6S 


6S 




^ 5/^(x) 


6S 


6S 


6r^{x) 


5l-{x) 



The right- and left-handed sources allow the definition of covariant derivatives 



D 



R 



ncolor _ • 



ncolor _ -7 



(A16) 

(A17) 
(A18) 

(A19) 
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where _D™'°'' is the covariant derivative with respect to local SU{3) color transformations. 
By construction, Dj^ip^i, Dj^ipi transform as the chiral fields themselves, and C is simply 
given by 

£ = i^^^^D^ijR + i^Lll^Dl^L + ^^Gauge , (A20) 

where /^cauge is the gauge field part of the QCD Lagrangian containing the gluon fields 
Obviously, flA20p is invariant under local chiral transformations. 
Having C in hand we can now define a generating functional for correlation functions 
involving the currents and densities, 

Z^cAr,.l,.sM = P[G^,^,^]e-/'^'^'^[^>'^-^'^], (A21) 

where Zq denotes the partition function, i.e. Zqqb for vanishing sources. Taking functional 
derivatives with respect to the sources one generates all possible correlation functions involv- 
ing the currents and densities; the generating functional for connected correlation functions 
can be defined as usual by the logarithm of Zqqd- 

The matching to the chiral effective theory is now done by requiring that the generating 
functional in the chiral effective theory is the same as the one in the fundamental theory, 

^QCD[r^,//.,S,p] = Zchir^i[rf,Jf,,S,p], (A22) 

Zctir4r^.Jf.,s,p] = ^ [ l?[vr]e-^'^*^^='^-'[^-'-^'Pl (A23) 

^0 J 



The equality in flA22p is meant in the sense that the two sides coincide order by order in 



a low energy expansion, where the long-distance correlation functions are dominated by 



the pion pole [29(. Correlation functions in the effective theory are obtained by the same 
functional derivatives as in the underlying theory. Provided the generating functionals, as 
functions of the sources, are the same, arbitrary correlation functions also agree. Moreover, 
the QCD symmetries are carried over to the effective theory, i.e. the Ward identities of 
QCD are correctly reproduced by the effective theory. Consequently, the right hand side 
must also be invariant under local chiral transformations, and this invariance provides one 
constraint in the construction of £chirai- Other constraints are provided by parity (P) and 
charge conjugation (C). 

The effective Lagrangian ^chiraii'^A*, ^n, s,p] is obtained in a systematic low energy expan- 
sion, and it has been derived by Gasser and Leutwyler through next-to-leading order 0] . 
We do not repeat their result here but emphasize that the local chiral invariance implies that 
the sources for the right and left handed currents can only enter the effective Lagrangian 
through the covariant derivative 

D^J: = d^T. + zSr^ - z/^S (A24) 
on the usual chiral field S, and through the field strength tensors 

r^v = d^Ty - dyT^ + z[r^, r^], 1^^ = d^ly - dj^ + . (A25) 
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Strictly speaking, the scalar source should not vanish but be set to the physical quark mass. 
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The transformation behavior of these objects is as expected, taking into account S L'LR\ 

L{D^T)R^ , r^, ^ Rr^,R^ , l^, ^ L/^.L^ . (A26) 

2. Generating functional in the Symanzik effective theory 

Following the development in the continuum, we want to define a generating functional 
in the Symanzik effective theory, which may then be matched to the chiral effective theory. 
As before we want to obtain the currents and densities by taking derivatives with respect 
to the source fields. As we will see this is not as straightforward as in continuum QCD. For 
simplicity we will deal here only with the currents and ignore the densities. This reveals the 
main obstacles in the procedure. 

As in the continuum case we want to define a source term in the Lagrangian, such that 
taking the derivatives with respect to sources produce the Symanzik currents, given in (fTTj) 
and f|T2|) . In addition, we want to do this in a way that maintains invariance under local 
chiral transformations, since this symmetry provides one of the links in the matching to the 
chiral effective theory. 

It is illustrative to first show that the naive generalization of the continuum case fails in 
this respect. Suppose we write down a source term that directly couples source fields to the 
vector and axial vector current, 

/^Source.Sym = (a^) V;^Sym,Loc (a^) + (^^i.^) ^Uju^^M ■ (A27) 

Just as in the continuum — cf. eq. (]A4p and (lASP — functional derivatives with respect to 
f°(x),a"(x) produce the desired currents. However, local invariance is lost: The 0(a) cor- 
rections in the currents and the derivative terms in (lAlOp . (lAlip . generate 0(a) terms under 



local chiral transformations that are not cancelled by the variation of the Symanzik effective 
action in ( fTOi) . That is of course no surprise, it simply reflects the fact that chiral symmetry 
is explicitly broken in the Symanzik effective theory and that the currents in ( fTTj) and ( |T2|) 
are not Noether currents associated with exact chiral symmetries. 

In order to maintain local chiral invariance we proceed differently and introduce separate 
sources fields for the 0(a) corrections in the currents. Our source Lagrangian has the form 

•^Source, Sym -^Source, ct ~l~ -^Source, a • (A28) 



The first term with subscript "ct" is the familiar one from continuum QCD, defined in (lA2p . 
As before, one can rewrite it in terms of right and left handed sources using flA15|) . The 
second part reads 



R 



+ {D'^i^R)l^C\lA,^ijL + ^Ri,rA,^C\D]:ijL . (A29) 
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This is essentially nothing but a source term for the 0(a) corrections to the currents as 
given in (!30l) and (!3T]) . We have introduced sources r^^, Ix fj_, X = V or A, which transform 
according to 



X,fi 



X,n 



I'. 



X,fi 



(A30) 



Note that this is the transformation law under local transformations even though there are 
no terms involving derivatives of L and R. Note also that i2sourcc a contains the covariant 



derivatives D^, D^, defined in ( ]A19l) . which involve the sources lij,,r^ for the continuum 



parts in the Symanzik currents. 

It is easily checked that the source term flA29p is invariant under local chiral transfor- 
mations (taking into account the transformation laws for Cv,Ca given in fl29l) ). Therefore, 
also the total Lagrangian 

(A31) 

is locally invariant. The currents are then obtained from S = J d^xC{x) according to 



5S 



6S 



6S 



+ 



6S 



5r-{x) 6l-ix) 5r-Jx) 5l'^^^{x) 



A 



/i,Sym 



5S 



Sr'^Ax) 



6S 
61^] 



+ 



6S 



6S 
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A,fi 



(A32) 



(A33) 



which is a generalization of the prescription in the continuum, (1A17P and (lAlSp . Note that 
we implicitly assume here that the spurion fields A, Cy, Ca are set to their "physical" values 
acsw, ocv, ac-A after the derivatives have been taken. 

So far we focused on the invariance under local chiral transformations. For the matching 
to the chiral effective theory the discrete symmetries P and C are also needed. In addition, 
there is one property in the source term that we also have to preserve. The source term in 
(]A29P depends on the sources and the spurion fields Cy and Ca- The dependence is such 
that the vector current sources ry, W only couple to Cy and the axial- vector sources only 
to Ca- a mixed term, involving C^^y^^ for example, is not present. It is mandatory to 
preserve this feature, otherwise the vector current in the effective theory may end up with a 
contribution proportional to ca, which is certainly not the case. In order to achieve this we 
choose to generalize the source term without changing the results for the currents derived 
from it. 

Instead of one spurion field Cy we introduce two of them, Cv,i and Cy^2-i and make the 
replacements 



Cyry, 



Cyiry 



W,iiCy 



W,iiCy^2 ' 



(A34) 



The symmetry properties of £source,a remains unchanged if both Cv,i and Cy^2 transform 
as Cy, and if both have the same "physical" value acy. However, the source term is now 
invariant under more general symmetry transformations, namely 

Hyry^fj_R\ Cyi > LCyiHy 

^ Lly,fj,Gy, 



I 



Cy2 



GyCyoR^ ■ 



(A35) 
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where Hy and Gy are two independent local SU{2) matrices. The origin of these two hidden 
local symmetries is the fact that the right and left-handed sources and the spurions Cv,i only 
appear next to each other in the source term. Note that we recover the previously discussed 
transformation laws with Hy = R and Gy = L. 

The same generalization can be done with the axial-vector part in the source term after 
we have introduced Ca,i and Ga,2 and postulated the transformation behavior 

fA,pL — * HArA,fiR\ Ca,i — > LCa,iH\ , 

U,/^ — ^ LlA,iiG\, Ca,2 — ^ GaCa.2R^ ■ (A36) 

Since the symmetries Ha-, Ga are independent of the ones in the vector current part we can 
no longer form invariants with vector current sources and axial-vector spurions Ca,! and 
vice-versa. This will become important in the next section where we construct the currents 
in the chiral effective theory. 

In analogy to flA21l) we can now define a generating functional in the Symanzik effective 
theory. The number of external fields, however, has grown significantly, 

^Sym = ^Sym['^/X5 ^^(5 ^^1,^(5 ^5 C*\/,25 C'a,!, C*A,2] 5 (A37) 

and here we still have not included any sources associated with the scalar and pseudoscalar 
density. Apparently, the method of constructing correlation functions via a generating func- 
tional looses its simplicity away from the continuum limit. In the naive continuum limit 
(1A37I) reduces to ^qcd with its dependence on r^, only. 



3. Matching to WChPT 



Similarly to the continuum case we now match the Symanzik effective theory to Wilson 
ChPT by demanding that the generating functionals in both theories agree 

•^Sym -^chiral ; 

(A38) 

where the right hand side is defined in analogy to ( 1A23|) . For simplicity we have dropped 
the dependence on the large number of external fields, specified in ( ]A37[) . 

The effective Lagrangian entering the right hand side in ( 1A38|) is build in terms of the 
pseudoscalar field S and its covariant derivative D^Tj as well as all external fields. Many 
terms have already been constructed and can be found in the literature. Firstly, one finds 
the familiar Gasser-Leutwyler Lagrangian [l|, [2I, which provides the continuum part. Terms 
involving the spurion field A are given in Ref. Hence we focus here on the new con- 
tributions stemming from the 0(a) terms in the currents, i.e. those involving the fields 

'^v,^l1 '^A,^i lA,fi, Cv,i, Ga,!- 

It will be useful to introduce the following combinations: 

Sv,^l = i{Gv,ifv,ii + lv,tiGv.2) ■, Sa,^i = i{GA,irA,fj. + lA,fj.GA,2) , 
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Fiplrl 


( ,h ^vs^ I 




L CXil 1 L y J 


















A 




/i 


/It 








I ^) ^X,ii 




Lrv g\. 


' X,ti 




Cx,i 


LCxH^x 


^X,l 


'-■x.i 


Cx,2 


GxCxR^ 


^X,2 


'~'X,2 


Sx,n 


LSx,,.R^ 


_qT 




Ax,,, 


LAx,,,B) 




(-1)^4,. 



TABLE I: Summary of transformation properties. The shorthand notation (—1)^ in the vector 
quantities represents the sign flip in the spatial components under parity. 



These combinations are automatically invariant under the four hidden symmetries involving 
HxiGx with X = V,A. The factor i is just convention, inspired by the observation that 
the terms involving the continuum sources also always include an i. For convenience we 
have summarized the transformation behavior of these quantities in table [H We also list 
the transformation rules under parity and charge conjugation, which one also needs for the 
construction of invariants under all symmetries. 

We are now in the position to construct the terms of 0(a) that will contribute to the 
currents. For this we need invariants involving one power of the sources in flA39l) . Lorentz 
invariance then requires one covariant derivative The simplest invariants that can be 
formed are 

(D^SA^^^ + AxAD^^^V) , X = V,A. (A40) 

Analogous invariants involving Sx,fi cannot be build since they violate charge conjugation. 

We do not need to consider any invariants involving D^Sx,,,, D^Ax^^i or terms quadratic 
in Sx,,i or Ax,,i- Even though there exist non-vanishing invariants they will not give con- 
tributions to the vector or axial- vector current. The reason is that all these terms contain 
more than one source field and derivatives of Cx- Setting in the end the sources to zero 
and Cx to its constant final value all these terms vanish. We also do not consider the terms 
involving (T) D^H ± S(D^S)"'') since 

tr (S^LI^S) = S^lndetS = 0. (A41) 



It turns out that the two terms given in flA40p are the only independent invariants at the 
order we are interested in. 

In the last section we stressed the importance of the local hidden symmetries for the 
correct chiral Lagrangian. If instead of (1A35I) . (1A36I) we assume the weaker transformation 
laws ( 1A30I) and ( !29l) . we can construct the additional invariants 

(StD^Srx,^ + S(Z}^S)t/x,^)(CxSt + ^C\) (A42) 
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These invariants require to "split" the right and left handed source fields from the 0(a) 
spurious Cx and build products of two trace terms, something that is forbidden by the 
transformation rule (]A35p . (1A36I) . 

It is straightforward to compute the 0(a) corrections to the currents following from (]A40p . 
In the case of the vector current the correction vanishes. For the axial-vector current we 
find the correction 

2acA9^(T"(S - St)) (A43) 

In order to derive the complete currents including all 0(a) corrections we also need other 
terms which have already been given by Ref. [sl. Carrying over their notation, the relevant 
terms in the Lagrangian are 

C = ^(D,,S(D^S)t) + 1^45(D,,S(D^S)t)(ASt + SAt) 

+W^o{D^nD^.A)^ + D^A{D^T.)^) + 2Wa{D^T.A\^ + A^,^(D^S)t) . (A44) 

The last term, coming with a new LEG Wa-, is the one we found above. We do not include 
the term with X = V since its contribution to the vector current vanishes anyway. 

Before deriving the currents we are free to make a field redefinition in order to simplify 
C. Following Ref. we perform the change S — >^ S + (5S with 



5S = -^(SAtS-A) (A45) 



and obtain C in (1A44P with the modified coefficient W45 W45 + Wiq. Therefore, without 



loss we can drop the Wiq term in (1A44P . as long as we consider only physical quantities. 

In fact, if we work only to linear order in the sources, it is easy to see that the Wa 
and WiQ terms in (1A44P are proportional, so that Wa can also be absorbed into W45 by a 
change of variables [9|. This holds as long as we consider only physical quantities and if all 
currents are placed at different space-time points. In order to avoid the latter restriction, and 
maintain generality, we do not make this change of variables. Nevertheless, for the quantities 
we calculate in this paper, which are both physical and involve separated currents, the 
possibility of this change of variables implies that the LECs must enter in the combination 
W^45Csw + WaCa- This provides a check on results. 

We can now compute the vector and axial- vector current according to (]A32[) . (]A33I) . and 



find exactly the same expressions given before in eqs. ( l43i) and (H4l) . We therefore conclude 
that we obtain identical results for the currents with the generating functional as with our 
"direct method". 



APPENDIX B: COMPUTATION OF Za loc- 

In this appendix we derive the result ( l60l) for the renormalization constant Za,loc, which 
follows from imposing the Ward-Takahashi identity (WTI) in (1591) . It will be useful to first 
show that fl59l) is indeed an identity in (massless) continuum ChPT. The generalization to 
WChPT is then straightforward. 
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We start by establishing 



at leading order in continuum CliPT. The operator Cout is composed of fields with support 
outside the time interval [yo — t,yo + 1]. We will not need to specify Oout in detail until we 
consider the 0(a) corrections; for now, we only assume that it creates an even number of 
pion fields, that the two-pion component occurs at LO in a chiral expansion, and that the 
two-pion component contains no zero-momentum pions. 

Expanding the vector current to leading order in the pion fields. 



ya 



fB2) 



(using abbreviation vr^ = S^tt"), the right-hand-side of the identity we want to show is simply 
given by 



rhsLo = 2i{VQ^^o{y)Oont)LO 

= -2e^'^\n^{y)4{y)Oout)LO 



(B3) 
(B4) 



where (C?)lo denotes functional integrals with -C^iji-ai = ^^/2 in the Boltzmann weight, and 
interactions, which stem from the expansion 



chiral 



1 



(B5) 



used at the lowest order giving a non- vanishing result. In fact, to evaluate rhsLo we need 
no interactions, given that Oout has a LO two-pion component. 

For the left-hand-side we need the expansion of the axial- vector current. 



^M,LO = y(r'^(St9,S-S9,St))=z/ 



0{iT^) .(B6) 



The contribution with the leading order term in both axial- vector currents on the left-hand- 
side vanishes, 



- f / ({<(2/o + t,x)- n^iyo -t,x)} n'o{y)Oout) 



LO 







(B7) 



since Oout does not contain zero momentum pion fields. 

The non-vanishing contribution to the left-hand-side stems from the three-pion term in 
one of the axial- vector currents or the four-pion term in £chirai- 

We obtain the first contribution using the three-pion term in the current Ao(y), finding 



Ihsi 



LO = -le^'" I dx{{7r-,{yo + t,x) 



KiVo - t,x)} {tt^tt ■ TToM - nyiy)} Oout) 



LO 



(B8) 



29 



Performing the Wick contractions this can be written as 

lhsi,Lo = 2e-''(7r«(y)7ro^(y)aut)LO J dx do {G{t, x - y) - Gi-t, x - y)} , (B9) 



where G denotes the massless pion propagator 

A = A"* / 

(27r)4p' 



(7r'^(x)7r^(y))Lo = S^'Gix-y) = 6'^' I ^^^M-v) . (BIO) 



The integral over x is easily evaluated. Using 

dxdoG{x) = [-} (Bll) 
I +^ tora;o < 

the integral reduces to —1 and we find that 

lhsi,Lo = rhsLo • (B12) 

We obtain the second contribution to the left hand side using the three-pion term for 
Ag(|/o ±t,x), yielding 



lhs2,LO = 2e'^''^y rff({7r"(y+)7r^(y+)9oG(t,f-f)-7r'^(y_)7ro^(y_)9oG(-t,f-y)}aut)LO, 

(B13) 

where we introduced the shorthand notation y± = [yo ±t,x). 

The third contribution stems from the four-pion terms in the chiral Lagrangian (1B5I) , and 
we find 

lhs3,LO = le^'"! dx {{7T^oiy+) - <iy~)} J rfS{(7r-7r^)2-7r2vrJ}(^) 7ro(y)aut)LO • 

(B14) 

This reduces to 

lhs3,Lo = e'^''^ j d^zdx (7r«(^)7rJ(^) {d^doG{y+ - z) - d,doG{y^ - z) 

+doG{y+ - z)d, - 9oG(y„ - z)d,} doG{z - y)0,^C)^o ■ (B15) 

Since 

j dxd^doG{y± - z) = -5^o5(z/o ±t- zq), (B16) 
the first two terms in (IBISP become 



-e'^'^ I c?x({7r'^(y+)7ro^(y+)9oG(t,x-y) -7r'^(y_)7r^(y_)9oG(-t,f -y)}aut)LO 

(B17) 
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where we renamed ^ as x. Using (IBlip for the integral over doG{y± — z) the remaining two 
terms become 



Jyo-t J 

and after partial integration this can be written as 

-e-'^ ! dx{{Tx\y+)Til{y+)d,G{t, x - y) - Ti'^{y^)nl{y^)d^G{-t, x - y)] O^^Oi^o 



ryo+t r 

+6'^^^ / dzo / dz{n-{z)d,d,Az)doG{z - y)aut)LO • 

Since df^7i'^{z) is contracted with the on-shell fields in Oout (recall our assumption!), the 
second line vanishes. The remaining first line is the same as in ( ]B17I) . Hence, in total the 
third contribution reduces to 

lhS3,LO = -lhS2,LO, (B18) 

so the second and third contributions cancel. This, together with eq. (1B12I) . proves the 
Ward-Takahashi identity fIBip to first non-trivial order in the chiral expansion. 

Repeating this calculation with the lattice spacing corrections included is now straight- 
forward. What changes are the expansions of the currents and the effective Lagrangian in 
terms of the pion fields: 

CcMr.i = ^(1 + X,)nl + ^ [(1 + X^){n ■ n^)^ - (1 + X2)7i\] + 0(71^) , (B19) 
V;% = ^e'^^^V^ll + Xi) + 0(7r^) , (B20) 
A;:,,^ = tf{l + Y,)nl + 2^ [(1 + Y,)it'^{t, ■ vr,) - (1 + Fs^vr^] , (B21) 

where we introduced the coefficients 

16 

Xi = -pH^45acsw, 
40 

X2 = j^W^^dcsw , 
g 

Fi = Xi + —WAdcA , 
Y2 = X.-^Wacica, 

Y3 = '^W^^dcsw + ^WAdcA . (B22) 

The modification of the kinetic term by the factor 1 -|- Xi yields a slightly different pion 
propagator. 



Sjab Sjab 
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with G defined in (IBlOp and (C)Lo(a) denotes functional integrals with C]:^^^^'^ = {l+Xi)ir^/2 

in the Boltzmann weight. Notice that 1 + Xi = [cf. (1581) ]. Hence, in terms of the 
renormalized pion fields 

rix) = Z-'/'n^x) , (B24) 
equation flB23|) assumes its standard form 

(7f«(x)7f^(y))Lo(a) = S'^'Gix - y) , (B25) 
The right hand side of the WTI now reads 

rhsLo(a) = -2e-''\l + (vr'^(y)vrJ(y)aut)LO(a) • (B26) 
while the three contributions on the left hand side are modified as follows: 

[1 + (3Yi + Y2 + 2^3) /3] o^cabw fe/„vo \ (V'07\ 

lhsi,LO(a) = {l^Xi) ^ (y)7r^(y)Cout)LO(a) , (B27) 

[i + (3yi + y2 + 2F3)/3] „ . 

lhs2,LO(a) = {l + Xi) ^ ~ lhs2,-) , (B28) 

lhs2,± = 2e"^^y" rfx9oG'(±t,x-y)(7r'^(y±)4(y±)aut)LO(a), (B29) 

lhS3,LO(a) = ^^^^^2 (^^^2,+ - lhS2,-) , (B30) 

Here the factors in the numerator come from (]B22p while those in the denominator arise 
from the presence of (one or two) pion propagators in the manipulations given earlier in the 
appendix. 

To determine the renormalization constant we need the ratio of the two sides of the 
original WTI. It follows from (lB26|l and (!B27l) that 

lhSi,LO(a) _ [1 + (3^1 + ^2 + 2^3)73] 

rhsLO(a) (l+^l)^ ' ^ ^ 

independent of the detailed form of Oout- 

For the remaining ratios we need to specify the external fields. We take 

OL, = e^''n\T,p)r{-T,-q) (B32) 

with p, q ^ 0. Here -7f'^(T, j?) is the Fourier transform of -7r'^(T, x) with respect to the three 
spatial coordinates. For this choice, the correlators we need are 

rhsLo(a) = -12 d^{G{T - yo,p)G{T + y^, q)} , (B33) 

lhs2,± = 12 ¥l^l^Zllo^{G{T -yoT t,P)G{T + yo±t, q)} , (B34) 
(i + Aij 

where = d/dT. For simplicity, and without loss of generality, we have set y = 0, which 
avoids an overall phase. 
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To evaluate these expressions we need the hybrid position-momentum pion propagator, 



using which we find 



EpEq 

lhs2± = ±1 3(Ep + Eg) ^^_\t\E^_^_^E^j^E„)T+{yo±t){E,-E„)] ^ ^^^^^ 

(1 + Xi) 2EpEq 

where Ep^g = |p — g|. Inserting these results into eqs. (]B26p . (1B28P and (]B30P we obtain 

— = :Y^2 cosh[t{\p\-\q\)]exp[-\t\\p-q]], (B38) 

lhs3,LO(a) _ [1 + 2^1 + {X, + 2X2) /3] 
rhsLo(a) (1 + ^1)^ 

and thus, finally. 



cosh[t(|p1 - \q\)]exp[-\t\\p-q\], (B39) 



IhSLO(a) _ lhSi,LO(a) + lhS2,LO(a) + lhS3,LO(a) 
rhSLO(a) rhsLo(a) 



(B40) 



= l + ^(iy45Cgw + ty^cA){l-cosh[t(|p1-|g-1)]exp[-|t||p-g-1]} . (B41) 

This ratio should be unity when multiplied by the renormalization factor (Za,loc/^a)^' 
leading to the result fl60l) in the main text. 
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